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Abstract 

By performing an analysis in moment space using high statistics DIS world data, we extract the values 
of both the QCD parameter A.^L up to NLO and of the power corrections to the proton structure 
function, F2 . At variance with previous analyses, the use of moments allows us to extend the kinematical 
range to larger values of x, where we find that power corrections are quantitatively more important. 
Our results are consistent with the n dependence predicted by IR renormalon calculations. We discuss 
preliminary results on nuclear targets with the intent of illustrating a possible strategy to disentangle 
power corrections ascribed to IR renormalons from the ones generated dynamically e.g. from rescattering 
in the final state. The latter appear to be modified in nuclear targets. 
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By performing an analysis in moment space using high statistics DIS world data, we extract the values of both 
the QCD parameter A^L up to NLO and of the power corrections to the proton structure function, F2. At 
variance with previous analyses, the use of moments allows us to extend the kinematical range to larger values of 
x, where we find that power corrections are quantitatively more important. Our results are consistent with the 
n dependence predicted by IR renormalon calculations. We discuss preliminary results on nuclear targets with 
the intent of illustrating a possible strategy to disentangle power corrections ascribed to IR renormalons from the 
ones generated dynamically e.g. from rescattering in the final state. The latter appear to be modified in nuclear 
targets. 



1. QCD FITS TO MOMENTS OF F 2 

The Q 2 dependence of the structure functions 
in Deep Inelastic Scattering (DIS) as well as in 
other related processes (e.g. Drell-Yan) can be 
described accurately within perturbative QCD 
(pQCD) up to NLO provided one confines calcu- 
lations to the kinematical region of Q 2 > 10 GeV 2 
and Bjorken x: 0.1 < x < 0.7 (correspond- 
ing to high values of the invariant mass squared: 
W 2 > 4 GeV 2 ). At the border of this kine- 
matical region the agreement between pQCD cal- 
culations and experiment is no longer accurate 
nor unambigous because of the increasing impor- 
tance of non-perturbative (np) physics. At large 
x the agreement with the data is known to im- 
prove by including power corrections 0J^] . Aside 
from parton model based phenomenological in- 
terpretations, the presence of power corrections 
still constitutes a rather elusive problem for the- 
ory. Canonical methods such as the hard scat- 
tering factorization (HSF) or OPE (for a review 
see H and references therein), are hampered by 
the large and hard-to-classify number of opera- 
tors. Comparison with data does not allow yet to 
distinguish among different higher twist contri- 
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butions (see e.g. (J]). More recently the idea has 
been developed that a one to one correspondence 
can be defined between the dynamically gener- 
ated higher twist terms and the ambiguities in the 
resummation of the asymptotic pQCD series, or 
renormalons (|| and references therein). As for 
np methods, there seems to be indications from 
lattice studies that power corrections do surface 
though not seemingly directly related to renor- 
malons Q. 

A preliminary step in order to test recent de- 
velopments is to explore whether, given the ac- 
curacy and kinematical coverage of present data, 
power corrections can be separated unambigously 
and their actual size determined, from NLO and 
higher order terms, and from a number of other 
non-QCD contributions. With the aim of achiev- 
ing a clear-cut answer, we started our analysis 
from the most well known structure function, F 2 . 

We constructed moments, 

M n (Q 2 )= f dxF 2 {x,Q 2 )x n - 2 , (1) 
Jo 

from the NMC, BCDMS and SLAC data sets. 
Differently from previous analyses, we set no 
lower limit on the final state invariant mass, thus 
extending the kinematical domain in x up to 
x ks 1 and Q 2 > 2 GeV 2 . Moments are more di- 
rectly connected to OPE, they can be calculated 
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Table 1 

Comparison of world measurements of Aqcd = A^l^_ and the coefficient t 2 



Collaboration 


Range in Q 2 GeV 2 


A^MeV 


r 2 Gel/ 2 


Before 1980 [3| 


1.9-85 


300 <-> 1000 


-0.1 <-> 0.25 


NMC § 


0.5-260 


250 


0.18 ±0.12 


NMC+BCDMS+SLAC [2] 


0.5-260 


263 ± 42 ± 55 


0.09 ±0.04 


CCFR [11] 


5 - 199.5 


371 ±31 


0.24 ±0.1 


BEBC-WA59 [12| 


Q 2 < 64 




-0.16 ±0.05 


This paper 7 4 = 


5 - 260 


241 ± 36 


0.21 ±0.09 


This paper 7 4 ^ 


2 - 260 


250 


0.25 ±0.20 



in principle using non-perturbative methods and, 
on the practical side, they are the only way to 
include data at large x without getting into the 
complications of the resonance structure. 

We performed a fit of the following expression 

M„(g 2 ) = Mp™ + £ A4 2fe + 2 >(g 2 ) -L , (2) 

where the Q 2 dependence of M$® GD = 
M n (p?)C n {Q 2 , [i 2 ), fi 2 being a given scale, was 
considered up to NLO. Furthermore, since we 
are interested in the large x behavior, we re- 
stricted the analysis to the Non Singlet (NS) con- 
tribution therefore avoiding ambiguities associ- 
ated with the gluon distribution. The subtraction 
of the Singlet contribution introduces some error. 
However, based on current parametrizations we 
estimated it to be 7—12% of the total moment for 
n = 3 and Q 2 = 2 — 50 GeV 2 and to be less then 
1% for n = 8. Further sources of systematic er- 
ror arise from the extrapolation to regions where 
no data are available which constitute 11.3% and 
23.5% of the whole kinematic domain contribut- 
ing to M n at n = 3 and n — 8, respectively, and 
from Target Mass Corrections (TMC) which dis- 
play an inverse power-like behavior and increase 
with n 0. In order to be able to compare with 
previous extractions we will show here results ob- 
tained in the "factorized" approximation, i.e. 

M n (Q 2 ) = M^ CD (l + + ai?^j , (3) 

corresponding to M„ = a„ t 2 etc., where the n 
dependence of the power corrections is included 
in the functions an , i — 1,2; t 2 and are the 



parameters to be determined by the fit, besides 
the pQCD ones; no Q 2 dependence is assumed 
for M n 4 ^ . Power corrections of order higher than 
0(1/ Q 4 ) are found to be negligible for Q 2 > 
2 GeV 2 . In order to compare with similar re- 
cent extractions pH, we set in Eqs.(^) and (0): 
7 4 = 0, and t 2 = M n 4) {Q 2 )/nMPQ CD (Q 2 ). Our 
results can be also anti-Mellin transformed to re- 
construct F 2 (x, Q 2 ), by using the properties of the 
Mellin transforms. Since the n dependence found 
in different models is rather simple, no polyno- 
mial technique was found to be necessary at this 
level. 

2. DISCUSSION OF RESULTS 

Our results for the parameters A-^gL and r 2 
are given in Table and compared with the 
ones obtained in previous extractions. Since re- 
cent analyses used directly F 2 (x,Q 2 ), giving IE- 
dependent power correction coefficients, Cht{x), 
we calculated the values of r 2 shown in the Ta- 
ble by performing Mellin transforms of Cht{x), 
times the correspondong pQCD structure func- 
tions parametrizations. The form a„ = n, trans- 
forming into Cht{ x ) ~ 1/(1 — x) at large x was 
used consistently. 

The first line of Table @ was included in order to 
show the highly enhanced accuracy provided by 
the NMC, BCDMS and SLAC set of data where 
it is now possible to perform a more careful ex- 
traction of power corrections. An anti-correlation 
clearly emerges from our fits between the param- 
eters of pQCD, e.g. A, and the coefficients of the 
power corrections. We studied this correlation; 
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details of our study are presented in || . Here it 
suffices to say that it originates from a region in 
Q 2 (5 < Q 2 < 20GeV 2 ) where neither the LO, 
simply logarithmic pQCD behavior, nor the pure 
power correction terms are dominant, whereas a 
mixture of higher order and 1/Q 2 terms are si- 
multaneously present. 

We dealt with this problem by extracting A us- 
ing only data at Q 2 > 30 GeV 2 and by including 
lower values of Q 2 stepwise in the fit, at the fixed 
value of A. The value of A we found is consistent 
with 250 MeV i.e. the value previously found in 
. The results presented in Table [l] are there- 
fore clear-cut in predicting that, for similar values 
of A, we find a contribution from power correc- 
tions which is about twice as large as in the pre- 
vious QCD analysis. We seem however to find 
a contradiction with the analysis using neutrino 
data in that these show both a larger value of 
A and a slightly enhanced value of t 2 . By in- 
cluding the 0(1/Q 4 ) term and by using our ap- 
proach of lowering the Q 2 threshold from above, 
we found out that j 4 — within error bars until 
one reaches the value of Q 2 nin = 5 GeV 2 . Finally, 
r 2 and 7 4 are also anti-correlated. By includ- 
ing in the fit both parameters, r 2 and 7 with 
the lowest allowed threshold, Q 2 min = 2 GeV 2 , we 
are able to fit the data with negative values of 
7 4 and r 2 0.25 GeV 2 , i.e. consistent with the 
value obtained at a higher threshold. We would 
like to point out however that due to the anti- 
correlation between the two parameters, by per- 
forming a general fit to the data regardless of Q 2 
"thresholds" , large positive values of 7 4 and small 
values of t 2 can be also found such results 
are to be interpreted as a manifestation of the 
anti-correlation between parameters. 

With a firm understanding of the different Q 2 
dependent contributions at intermediate and low 
Q 2 , one can subsequently investigate the renor- 
malon hypothesis. We present here results of our 
study of the n-dependence of the coefficient of 
the 1/Q 2 term, aiPr 2 , in Eq.(|). In Figure 1 
we repeated our fits using the forms: aj, = n, 
an** — n 2 , corresponding to the dominance of 
two-quark and four-quark higher twist diagrams, 
respectively, in e.g. HSF models and we com- 
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Figure 1. The value of r 2 vs. n extracted ac- 
cording to different models for the n-dependent 
coefficient a n (l) as explained in the text. 



pared them to fits using the n-dependence pre- 
dicted by renormalon calculations (for details see 
references [9] and [10] in ||). Since in all predic- 
tions the n-dependent term is factored out from 
an unknown constant term, one can discriminate 
"models" based on the "flatness" of the n behav- 
ior extracted from the data. From Figure 1 one 
can see that the renormalon model gives clearly 
a better interpretation than the widely used as- 
sumption, = n. 

3. CONCLUSIONS AND FUTURE DE- 
VELOPMENTS 

Our initial study of power corrections in DIS is 
meant to be an exploratory one centered around 
two points: we first addressed the question of 
whether there is at all a direct and clear evidence 
of power corrections in the data, and if their ac- 
tual magnitude can be extracted; we then moved 
on to a phenomeonological study of their nature, 
namely whether they can be associated to long 
distance features of pQCD (renormalon models) 
or if other dynamical mechanisms pertinent to 
DIS only, such as final state rescatterings should 
be taken into account. Our main conclusions on 
the first part are that indeed power corrections 
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can be separated out from present data by using 
an approach in which different portions of the Q 2 
domain (proceeding from high to low Q 2 values) 
are analyzed separetely. The magnitude of the 
0(1/Q 2 ) coefficient is r 2 « 0.2 GeV 2 , which will 
need to be systematically compared with power 
corrections extracted from other reactions in fur- 
ther investigations, to test the presence of a uni- 
versal scale. An important outcome of our study 
is also that, within the present accuracy and kine- 
matical coverage of the data, meaningful results 
should be presented including the correlation be- 
tween the different parameters of the fit. Our ex- 
tractions seem to be consistent with renormalon 
calculations. 

More detailed knowledge on the nature of 
power corrections can be sought in principle by 
using nuclear targets. By assuming the validity 
of IA, i.e. by disregarding nucleon rescattering 
effects, one obtains a factorization of the A and 
Q 2 dependences of the moments urn 



M*{Q*)=M? l M^{Q 2 ), 



(4) 



Deviations from such a factorization would cor- 
respond to to a breakdown of IA and they would 
therefore signal the presence of rescattering ef- 
fects. As shown in uM present data on nuclear 
targets are not accurate enough to observe any 
feature in the ratio of moments and they sup- 
port the IA hypothesis. More detailed studies, 
together with higher accuracy nuclear DIS data 
might be able to reveal an underlying structure 
0. At present, by adopting the factorized form, 
Eq.(^), we can extract the value of t 2 for a nu- 
cleon. Results are presented in Figure 2 and com- 
pared with the proton values. From the Figure 
we can see that the proton and nucleon results 
are consitent with one another, i.e. no difference 
is seen in the neutron contribution. 

I thank the Institute of Nuclear and Particle 
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